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Figure 1: Left panel: The phase diagram of QCD on the hypothetical light quark mass versus strange quark 
mass plane. Thick lines correspond to second order phase transitions, the purple regions represent first order 
phase transitions and the yellow region represents an analytic cross-over Right panel: The most popular 
scenario for the /i-T phase diagram. At low temperatures and low chemical potentials we have the hadronic 
phase, which is separated by a cross-over from the low density high temperature quark-gluon plasma phase. 
The cross over ends in a critical endpoint (E) after which a first order transition separates the different phases. 
The low temperature high density region is conjectured to be a colour superconducting phase. 

1. Introduction 

The QCD transition at non-vanishing temperatures {T) and/or baryonic chemical potential (/i) 
plays an important role in the physics of the early Universe and of heavy ion collisions (most 
recently at RHIC at BNL; LHC at CERN and FAIR at GSI will be the next generation of ma- 
chines). The main goal of the present summary is to present some selected results of the Budapest- 
Wuppertal group on the QCD transition at vanishing and non-vanishing chemical potential. Most 
of the T=0 results were obtained at four different sets of lattice spacings and a careful continuum 
extrapolation was performed, we consider them as full results. 

The standard picture for the QCD phase diagram on the light quark mass {mud) vs. strange 
quark mass (m^) plane is shown by the left panel of Figure [l|. It contains two regions at small and at 
large masses, for which the T > QCD transition is of first order. Between them one finds a cross- 
over region, for which the QCD transition is an analytic one. The first order transition regions and 
the cross-over region are separated by lines, which correspond to second order phase transitions. 

When we analyze the nature and/or the absolute scale of the T > QCD transition for the 
physically relevant case two ingredients are quite important. 

First of all, one should use physical quark masses. As the left panel of Figure |l| shows the 
nature of the transition depends on the quark mass, thus for small or large quark masses it is a 
first order phase transition, whereas for intermediate quark masses it is an analytic cross over. 
Though in the chirally broken phase chiral perturbation theory provides a controlled technique to 
gain information for the quark mass dependence, it can not be applied for the T > QCD transition 
(which deals with the restoration of the chiral symmetry). In principle, the behaviour of different 
quantities in the critical region (in the vicinity of the second order phase transition line) might give 
some guidance. However, a priori it is not known how large this region is. Thus, the only consistent 
way to eliminate uncertainties related to non-physical quark masses is to use physical quark masses 
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Figure 2: The ratio f{cont)lf{Nt) as a function of 1/N} (left panel). f{cont) is the continuum extrap- 
olated free energy of the staggered fermionic gas in the non-interactive, infinitely high temperature limit. 
f{Nt ) is the value obtained on a lattice with Nt temporal extension. The black line shows our choice (stout 
improvement, only next-neighbour terms in the action), whereas the red and blue lines represent the Naik 
and p4 actions, respectively. Masses and taste symmetry violation for different approaches in the literature 
(right panel). The smallest, physical quark mass and the smallest taste symmetry violation was reached by 
our works (black dot, |^ ^). Somewhat larger taste symmetry violation and about three times larger quark 
masses were reached by the MILC analysis on QCD thermodynamics (red dot, |[Tc|]). Even larger taste 
symmetry violation and about four times the physical quark masses are the characteristics of the Bielefeld- 
Brookhaven-Columbia-Riken result on Tc (blue dot, |Q). 

(which is, of course, quite CPU demanding). 

Secondly, the nature of the T > QCD transition is known to suffer from discretization er- 
rors [|l|, ^, |3|]. The three flavour theory with standard action on Nt=4 lattices predicts a critical 
pseudoscalar mass of about 300 MeV. This point separates the first order and cross-over regions of 
Figure [l| (left panel). If we took another discretization, with another discretization error, e.g. the 
p4 action and Nt=4, the critical pseudoscalar mass turns out to be around 70 MeV (similar effect 
is observed if one used stout smearing improvement and/or Nt=6). Since the physical pseudoscalar 
mass (135 MeV) is just between these two values, the discretization errors in the first case would 
lead to a first order transition, whereas in the second case to a cross-over. The only way to resolve 
this inconclusive situation is to carry out a careful continuum limit analysis. 

Since the nature of the transition influences the absolute scale (Tc) of the transition -its value, 
mass dependence, uniqueness etc.- the above comments are valid for the determination of T^, too. 

Thus, we have to answer the question: what happens for physical quark masses, in the con- 
tinuum, at what r^? To get a reliable answer we used physical quark masses on Nt=4,6,8 and 10 
lattices, which correspond to approximately 0.3, 0.2, 0.15 and 0.12 fm lattice spacings, respectively. 

It was conjectured that the physical point, for which the quark masses are tuned to their phys- 
ical value, is in the yellow, cross-over region. (One tunes the quark masses to their physical value 
by tunig the pseudoscalar masses -pion, kaon- to their physical value.) As we will see this conjec- 
ture turned out to be true. We show that the continuum extrapolated lattice result with staggered 
fermions is indeed a cross-over for physical masses. The existence of a cross-over transition at 
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Figure 3: The volume dependence of the susceptibility peaks for pure SU(3) gauge theory (Polyakov-loop 
susceptibility, left panel) and for full QCD (chiral susceptibility on Nt=A and 6 lattices, middle and right 
panels, respectively). 

vanishing baryonic chemical potential is one important necessaiy condition for the most popular 
/i-T phase diagram scenario (c.f. right panel of Figure |]). In this picture the cross-over region on 
the jX-T plane ends in a critical endpoint, after which a first order phase transition appears. 

The only systematic way to get quantitative informations about the above features of the phase 
diagram is lattice QCD, which is extremely difficult for non-vanishing chemical potentials. As we 
have seen it is of crucial importance to extrapolate to the continuum limit in a controlled manner. 

In the presentation [Q] published results of the Bielefeld-Brookhaven-Columbia-Riken Col- 
laboration ^ from Nt=A and 6 were shown (and some unpublished figures for A^t=8, which were 
obtained within the HotQCD Collaboration). Since the CPU requirements for thermodynamics in- 
crease as K, Nf^ our Nt=\0 simulations need about 50 times more CPU than Nt=(). Do we have 50 
times more resources for QCD thermodynamics than our competitors? Of course not (it is almost 
the other way around). Instead, reaching A^,=10 is a fine balance. It is partly related to the choice of 
our action (which will be discussed in the next section), partly to the arrangements of the financial 
resources. For instance, as A'^; increases, one needs more and more statistics. Thus the thermal- 
ization can be done only once on a relatively expensive, scalable machine, such as Blue-Gene/L, 
whereas a large fraction of the non-vanishing T simulations can be done on more cost effective 
devices such as personal computer graphics cards. A 2 years old model can accommodate N,=6 
lattices, on a one-and-a-half year old model you can put Nt=?> lattices and the one year old model 
can work with quite large A'f=10 lattices. It costs a few hundred dollars and can provide upto 30- 
60 Gflops sustained QCD performance. They are not easy to code, adding two numbers needs 3 
pages, but recently more efficient programming environments were introduced. Clearly, this type 
of hardware provides a very advantageous price-performance ratio for lattice QCD. 

2. The choice of the action 

The first step is to choose an action, which respects all the needs of a thermodynamic analysis. 
T=0 simulations are needed to set the scale and for renormalization. T>0 simulations are needed 
to map the behaviour of the system at non-vanishing temperatures. The action should maintain a 
balance between these two needs, leading to approximately the same uncertainties for both sectors 
(otherwise a large fraction of the CPU-power is used just for "over-killing" one of the two sectors). 
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Figure 4: Continuum extrapolated susceptibilities T"^ / inP' as a function of \I{T^V). For true phase 
transitions the infinite volume extrapolation should be consistent with zero, whereas for an analytic crossover 
the infinite volume extrapolation gives a non-vanishing value. The continuum-extrapolated susceptibilities 
show no phase-transition-like volume dependence, though the volume changes by a factor of five. The V-^oo 
extrapolated value is 22(2) which is 11(7 away from zero. For illustration, we fit the expected asymptotic 
behaviour for first-order and 0(4) (second order) phase transitions shown by dotted and dashed lines, which 
results in chance probabilities of 10^'^ (7 x 10^'^), respectively. 

We used Symanzik improved gauge and stout improved 2-1-1 flavour staggered fermions [^] (due 
to the stout improvement we have only next-neighbour terms in the fermionic part of the action). 
The simulations were done along the line of constant physics. The parameters were tuned with 
a quite high precision, thus at all lattice spacings the niK/ fk and niK/m^ ratios were set to their 
experimental values with an accuracy better than 2%. 

The choice of the action has advantages and disadvantages. As we will see the advantages are 
probably more important than the disadvantages. The left panel of Figure ^ shows the continuum 
free energy divided by its value at a given Nf. A related plot is usually shown by the Bielefeld- 
Brookhaven-Columbia-Riken collaboration as a function of Nt. Since in staggered QCD most 
lattice corrections scale with a^, which is proportional to l/N^, it is instructive to show this ratio as a 
function of l/N^, for our action, for Naik and for p4. Extrapolations from 4 and 6 always overshoot 
or undershoot. Clearly, the Naik and p4 actions reach the continuum value much faster than our 
choice, but the scaling appears quite early even for actions with next-neighbour interactions. 
Extrapolations from Nt and Nt+2 with our action are approximately as good as Nt with the p4 
action (which was tailored to be optimal for this quantity, namely for the free energy at infinitely 
large temperatures). In practice, it means that our choice with Nt=8,lO gives approximately 2% 
error for the free energy. In a balanced analysis you do not need more, because the corresponding 
lattice spacings 0.15 and 0.12 fm are most probably not fine enough to set the scale unambiguously 
with the same accuracy. (E.g. the asqtad action at N^IO, which corresponds about a=0. 12 fm 
lattice spacing, still has ?alO% scale difference between ri & /k-) Since the p4 action is almost 20 
times more expensive than the stout action, it is not worth to pay this price and improve one part of 
the calculation, which hinders you to reach a reasonable accuracy in another part of the calculation. 

This is the balance one should remember in thermodynamics. Indeed, taste symmetry violation 
should be suppressed for many reasons (setting the scale at 7 = 0, restoring chiral symmetry at 
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r > etc.) As it was argued above it is more important to improve on this sector of the calculation 
than on the infinitely high temperature behaviour. The left panel of Figure § shows the splitting 
between the Goldstone and the first non-Golds tone pions for the Bielefeld-Brookhaven- Columbia- 
Riken Collaboration (which is far beyond their kaon mass), for MILC and for us. Our small splitting 
is partly related to the stout improvement and partly to the cost issues, since smaller lattice spacings 
could be used, which resulted in smaller splitting. In addition, the cost issues allowed us to use in 
our finite T simulations physical quark masses instead of much larger masses. 

3. The nature of the QCD transition 

The next topic to be discussed is the nature of 
the QCD transition. Physical quark masses were 
used and a continuum extrapolation was carried out 
by using four different lattice spacings. The details 
of the calculations can be found in [j8|]. In order 
to determine the nature of the transition one should 
apply finite size scaling techniques for the chiral 
susceptibility % = {T /V) ■ (d^logZ/dml^,). This 
quantity shows a pronounced peak as a function of 
the temperature. For a first order phase transition, 
such as in the pure gauge theory, the peak of the 
analogous Polyakov susceptibility gets more and more singular as we increase the volume (V). 
The width scales with lA'^ the height scales with volume (see left panel of Figure A second 
order transition shows a similar singular behaviour with critical indices. For an analytic transition 
(what we call a cross-over) the peak width and height saturates to a constant value. That is what 
we observe in full QCD on Nt=4 and 6 lattices (middle and right panels of Figure We see an 
order of magnitude difference between the volumes, but a volume independent scaling. It is a clear 
indication for a cross-over. These results were obtained with physical quark masses for two sets of 
lattice spacings. Note, however, that for a final conclusion the important question remains: do we 
get the same volume independent scaling in the continuum; or we have the unlucky case what we 
had in the Introduction for 3 flavour QCD (namely the discretization errors changed the nature of 
the transition for the physical pseudoscalar mass case)? 

We carried out a finite size scaling analyses with the continuum extrapolated height of the 
renormalized susceptibility. The renormalization of the chiral susceptibility can be done by taking 
the second derivative of the free energy density (/) with respect to the renormalized mass (m^). 

We apply the usual definition: //r^ =-N^-[logZ{N„Nt)/{NtN^) -logZ{NM,Nto)/{NtoN^o)]- This 
quantity has a correct continuum limit. The subtraction term is obtained at T=0, for which simula- 
tions are carried out on lattices with N^o, Nto spatial and temporal extensions (otherwise at the same 
parameters of the action). The bare light quark mass (/«„(/) is related to by the mass renormaliza- 
tion constant mr=Zm-mud- Note that Z,„ falls out of the combination m^d^ldm^j.=m^^^d^ldm^^. Thus, 
m^j [Xi^s^Nt) — xi^sOjNto)] also has a continuum limit (for its maximum values for different Nt, 
and in the continuum limit we use the shorthand notation m^Ax)- 
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Figure 6: Melting curves of different natural fats. 



In order to carry out the finite volume scaling in the continuum limit we took three different 
physical volumes (see Figure 0). The inverses of the volumes are shown in units of Tc. For these 
3 physical volumes we calculated the dimensionless combination T'^/ni^A;^ at 4 different lattice 
spacings: 0.3 fm was always off, otherwise the continuum extrapolations could be carried out, 
which are shown on Figure 0. Our result is consistent with an approximately constant behaviour, 
despite the fact that we had a factor of 5 difference in the volume. The chance probabilities, that 
statistical fluctuations changed the dominant behaviour of the volume dependence are negligible. 
As a conclusion we can say that the staggered QCD transition at mu=0 is a cross-over. (Note, that 
an analytic, cross-over like transition appears in other sector of the standard model, namely for the 
electroweak transition see e.g. [11] and references therein). 



4. The transition temperature 

An analytic cross-over, like the QCD transition has no unique T^. A particularly nice example 
for that is the water-vapor transition (c.f. Figure ||). Up to about 650 K the transition is a first 
order one, which ends at a second order critical point. For a first or second order phase transition 
the different observables (such as density or heat capacity) have their singularity (a jump or an 
infinitly high peak) at the same pressure. However, at even higher temperatures the transition is 
an analytic cross-over, for which the most singular points are different. The blue curve shows the 
peak of the heat capacity and the red one the inflection point of the density. Clearly, these transition 
temperatures are different, which is a characteristic feature of an analytic transition (cross-over). 

There is another -even more often experienced- example for broad transitions, namely the 
melting of butter. As we know the melting of ice shows a singular behavior. The transition is of 
first order, there is only one value of the temperature at which the whole transition takes place at 
O^C (for 1 atm. pressure). Melting of butter^ shows analytic behaviour. The transition is a broad 
one, it is a cross-over (c.f. Figure ^ for the melting curves of different natural fats). 

' Natural fats are mixed triglycerids of fatty acids from C4 to C24, (saturated or unsaturated of even carbon numbers). 
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Figure 7: Temperature dependence of the renormalized chiral susceptibility (m^AX\jf\i//T'^), the strange 
quark number susceptibihty (Xs/T^) and the renormahzed Polyakov-loop (Pr) in the transition region. The 
different symbols show the results for A', =4,6,8 and 10 lattice spacings (filled and empty boxes for Nt =4- 
and 6, filled and open circles for Ni and 10). The vertical bands indicate the corresponding transition 
temperatures and their uncertainties coming from the T^^O analyses. This error is given by the number in the 
first parenthesis, whereas the error of the overall scale determination is indicated by the number in the second 
parenthesis. The orange bands show our continuum limit estimates for the three renormalized quantities as 
a function of the temperature with their uncertainties. 
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Since we have an analytic cross-over also in QCD, we expect very similar temperature depen- 
dence for the quantities relevant in QCD (e.g. chiral condensate, strange quark susceptibility or 
Polyakov loop). 

In QCD we will study the chiral and the quark number susceptibilities and the Polyakov loop. 
Usually they give different Tc values, but there is nothing wrong with it. As it was illustrated 
by the water-vapor transition it is a physical ambiguity, related to the analytic behaviour of the 
transition. There is another, non-physical, ambiguity. If we used different observables (such as the 
string tension, ro, the rho mass or the kaon decay constant), particularly at large lattice spacings 
we obtain different overall scales. They lead to different Tc values. This ambiguity disappeares in 
the continuum limit. According to our experiences, at finite lattice spacings, the best choice is the 
kaon decay constant f^. It is known experimentally (in contrast to the string tension or ro), thus no 
intermediate calculation with unknown systematics is involved. Furthermore, it can be measured 
on the lattice quite precisely. 

Figure shows the results for the chiral susceptibility, for the quark number susceptibility and 
for the Polyakov loop. Red, blue, green, and purple indicate Nt=^,6,^ and 10 lattices. Nt^A is 
always off, the rest scales nicely. The shaded regions indicate the continuum estimates. There is 
a surprising several sigma effect. The remnant of the chiral transition happens at a quite different 
temperature than that of the deconfining transition. It is quite a robust statement, since the Polyakov 
transition region is quite off the ;^-peak, and the ;^-peak is quite far from the inflection point of the 
Polyakov loop. This quite large differece is also related to the fact that the transition is fairly broad. 
The widths are around 30-40 MeV. 

Due to the broadness of the transition the 
normalization prescription changes Tc, too. It is 
easy to imagine why, just multiply a Gaussian 
by and the peak is shifted. That means using 
X/T'^ gives about 10 MeV higher Tc than our 
definition, for which a normalization was ap- 
plied. (Note, that for the unrenormalized % nT^ 
normalization is natural, whereas for the renor- 
malized % the natural normalization is done by 
r^. This kind of naturalness manifests itself as 
possibly small errors of the observable.) 

Figure ^ shows the difference between the 
Tc values obtained by the Polyakov loop and by the chiral condensate as a function of the lattice 
spacing squared. The blue band indicates the difference for the chiral susceptibility peak position 
for the T^ and normalization. Thus using the T^ normalization no difference can be seen for 
M=4 and 6, a slight difference is observed for M=8 and a reliable continuum extrapolation needs 
A^,=6,8 and 10. 




[fm^] 

Figure 8: Difference between the the Tc values ob 
tained by the Polyakov loop and by the chiral con 
densate as a function of a^. 



Our resuh on Tc and that of the MILC Collaboration (rc=169(12)(4) MeV Ill0|]) are consistent 
within the (quite sizable) errorbars. 

However, our result contradicts the recent Bielefeld-Brookhaven-Columbia-Riken result 
which obtained 192(7)(4) MeV from both the chiral susceptibility and Polyakov loop. This value 
is about 40 MeV larger than our result for the chiral susceptibility (for the Polyakov loop suscep- 
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Figure 9: Possible contributions to the 40 MeV difference between the results of Refs. [ph and 
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Figure 10: Continuum extrapolations based on Nt=A- and 6 (left panel: inconsistent continuum limit) and 
using Nt=6,S and 10 (right panel: consistent continuum limit). 



tibility the results agree within about 1.5 a). What are the differences between their analyses and 
ours, and how do they contribute to the 40 MeV discrepancy? The most important contributions 
to the discrepancy are shown by Figure ^ The first difference is, that in [||] no renormalization 
was carried out, instead they used the unrenormalized quantity x/^'^- Due to the broadness of the 
distribution this observable leads to about 10 MeV larger than our definition. The overall errors 
can be responsible for another 10 MeV. The origin of the remaining 20 MeV is somewhat more 
complicated. One possible explanation can be summarized as follows. In Ref. [||] only Nt=4 and 
6 were used, which correspond to lattice spacings a=0.3 and 0.2 fm, or a^^=700MeV and IGeV. 
These lattices are quite coarse and it seems to be obvious, that no unambiguous scale can be deter- 
mined for these lattice spacings. The overall scale in Ref. [Q] was set by ro and no cross-check was 
done by any other quantity independent of the static potential (e.g. fk). This choice might lead to 



an ambiguity for the transition temperature, which is illustrated for our data on Figure [10|. Using 
only Nt=4 and 6 the continuum extrapolated transition temperatures are quite different if one took 
ro or fa to determine the overall scale. This inconsistency indicates, that these lattice spacing are 
not yet in the scaling region (similar ambiguity is obtained by using the p4 action of |^). Having 
Nt=4,6,8 and 10 results this ambiguity disappears (as usual Nt=4 is off), these lattice spacings are 
already in the scaling region (at least within our accuracy). This phenomenon is not surprising 
at all. As it was already mentioned e.g. the asqtad action at Nt~lO (which corresponds to about 
a=0.12 fm lattice spacing) has s:::10% scale difference predicted by ri or f^. 

The ambiguity related to the inconsistent continuum limit is unphysical, and it is resolved as we 
approach the continuum limit (c.f. Figure The differences between the values for different 
observables are physical, it is a consequence of the cross-over nature of the QCD transition. 
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5. Equation of state 



In the case of the two previous sections (nature of the transition and Tc) the full results are 
available in the staggered formalism. This is ensured by the fact that results from three lattice 
spacings out of four were found to be in the scaling region (c.f. two points are always on a line, 
thus the real behaviour can be read off only by having at least three points). This allowed 
controlled continuum extrapolations. For the equation of state the situation is not yet satisfactory, 
our group have results only for two lattice spacings (see also the contribution of F. Karsch |Q] for 
results on the equation of state; they also used two lattice spacings with the p4 action). 

Our goal is to determine the temperature dependence of the pressure (or energy density, en- 
tropy, speed of sound etc.) of QCD with physical pion mass of 135 MeV along the line of constant 
physics (LCP). For large homogenous systems the pressure is proportional to the logarithm of the 
partition function (poclog[Z]). All other quantities can be determined by using the pressure. Since 
the partition function itself is difficult to determine the usual technique is to calculate its derivatives 
and integrate. Thus, for the normalized pressure one obtains 
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Figure 1 1 shows the result for the normalized 
pressure as a function of T/r^. The two curves 
represent Nt=4 and Nt=6. Clearly, these two sets 
of lattice spacings are not enough for a controlled 
continuum extrapolation. Usually three points in 
the scaling region is needed to ensure that. Thus, 
in a fortunate case one could obtain the full result 
with Nt=4,6 and 8, whereas a somewhat less for- 
tunate case would mean Nt=6,i and 10 (as it was 
needed for the transition temperature, too). Note, 
that A',=8 needs about 30 times more CPU power 
than Nt=6 and A',=10 needs an additional factor of 
about 15. 

An old and serious problem for QCD thermo- 
dynamics is the link between perturbation theory and lattice QCD. While available lattice results 
for the equation of state (both for pure gauge theory and for full QCD) end at around 5 • , standard 
perturbation theory converges only at extremely high temperatures (at 5 • Tc, the different pertur- 
bative orders can not even tell the sign of the deviation from the Stefan-Boltzmann limit). Until 
recently no link between the two most systematic methods of QCD, namely perturbation theory 
and lattice QCD, existed for bulk thermodynamical quantities. It is of extreme importance to close 



Figure 11: The normalized pressure as a func 
tion of T /Tc. 



the gap between these results (for a recent result on that see [ 12]) 



6. Results at non- vanishing baryonic chemical potential 

In the previous three sections results were presented, some of which can be considered as full 
ones (the nature of the transition and Tc). For these observables three lattice spacings in the scal- 
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ing region were used and controlled continuum extrapolations were possible. For another quantity 
(equation of state) only two lattice spacings were applied, therefore no controlled continuum ex- 
trapolated results are available. The situation is even worse for non-vanishing baryonic chemical 
potentials (/i). In this case traditional Monte-Carlo simulations do not work. The determinant of 
the Dirac operator turns out to be complex, spoiling any method based on importance sampling. 
Five years ago new techniques appeared (starting with the multi-parameter reweighting technique 
[HI 0]), which can give predictions also at /^/O. Nevertheless, QCD at non- vanishing chemical 
potentials is much more CPU demanding than the }X=0 case. Thus, continuum extrapolated results 
will be most probably not available in the next one or two years for most observables. (There is at 
least one important exception. Using the present methods and technologies one can determine the 
continuum curvature for the T-}X phase diagram at /x=0.) In this section a few illustrative examples 
of our jX^O results are shown. 

According to the most popular scenario a critical point is expected on the temperature versus 
baryonic chemical potential plane (for an alternative possibility see Ref. [15] and the talk of O. 
Philipsen in these proceedings [16]). Using the multi-parameter reweighting technique [13, [I^ we 
studied dynamical QCD with rif =2+1 staggered quarks of physical masses on Nt=4 lattices []T7|]. 
The results are summarized on Figure 12 . At vanishing chemical potentials there is an analytic, 
cross-over like transition between the hadronic and quark-gluon plasma phases. As we increase 
H the transition temperature decreases. At a baryonic chemical potential around 360 MeV the 
cross-over region ends and a second order phase transition is observed. For even larger chemical 
potentials the transition is of first order. In order to tell the difference between an analytic cross- 
over and a singular transition (first or second order phase transition) a careful finite volume analysis 
is needed, similar to that of Section 3. This result is quite a different picture than that of Ref. [jl^]. 
Note, however that the two results are not at all in contradiction with each other. On the level 
they can be directly compared both observe no strengthening of the transition in leading order of 
/I. The strengthening in our result, which is needed for the critical endpoint, appears in much 
higher orders. (Note, that the same multi-parameter reweighting can be used e.g. to determine the 
equation of state for non-vanishing baryonic chemical potentials [jT8|, p^].) 

Another interesting feature of the jX-T phase diagram is the conjectured colour superconduct- 
ing phase at large jj. and small T. The presently available techniques are working well for small to 
moderate chemical potentials slightly below and all the way above the transition temperature. In 
order to extend the applicability range, we used [20] the density of state method. «/ = 4 staggered 
QCD was studied on 6^ and 8 • 6^ lattices. The method is expensive, thus only small lattices were 
used. Nevertheless, these small lattices show an indication for a triple-point connecting three dif- 
ferent phases on the phase diagram. The triple point is around Hq^^ 300 MeV and T 135 MeV 
(the quark chemical potential Hq is one third of the baryon chemical potential). Note, that the tem- 
perature will move most probably to smaller values in the continuum limit, since for Nt=6 lattices 
the smallest T is 73 MeV (if nip fixes the scale). The mass dependence was checked, at small T the 
position of the transition did not depend on the pion mass. 

It is important to emphasize again that results at non-vanishing chemical potentials are ob- 
tained on coarse lattices and the most important uncertainty is the systematic error due to large 
lattice spacings. 
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Figure 12: The phase diagram on the temperature versus baryonic chemical potential plane. Left panel: 
locating the critical point on Nt=A- lattices. Right panel: application of the density of state method to separate 
the three phases in n/=4 QCD. In both cases the continuum extrapolation is missing. 



7. Conclusion 

Some QCD thermodynamics results from the Budapest-Wuppertal group were summarized at 
vanishing and at non- vanishing chemical potentials. The necessary balance between T = and T > 
simulations was discussed in detail. As a consequence, Symanzik improved gauge and stout-link 
improved staggered fermionic lattice action was used in the simulations with an exact simulation 
algorithm. Physical masses were taken both for the Ught quarks and for the strange quark. The 
parameters were tuned with a quite high precision, thus at all lattice spacings the mx / fk and mk / m^i 
ratios were set to their experimental values with an accuracy better than 2%. Up to four sets of 
lattice spacings on lattices with A^,=4,6,8 and 10 temporal extensions were used (they correspond 
to lattice spacing a~0.3, 0.2, 0.15 and 0.12 fm) to carry out the continuum extrapolation. It turned 
out that only Nt=6,8 and 10 can be used for a controlled extrapolation, Nt=4 is out of the scaling 
region. 

The nature of the T>0 transition was determined. The renormalized chiral susceptibility was 
extrapolated to vanishing lattice spacing for three physical volumes, the smallest and largest of 
which differ by a factor of five. This ensures that a true transition should result in a dramatic in- 
crease of the susceptibilities. No such behaviour is observed: the finite-size scaling analysis showed 
that the finite-temperature QCD transition in the hot early Universe was not a real phase transition, 
but an analytic crossover (involving a rapid change, as opposed to a jump, as the temperature var- 
ied). As such, it will be difficult to find experimental evidence of this transition from astronomical 
observations. Since for present day heavy ion experiments the baryonic chemical potential is also 
very small, the above results apply for them, too. 

The absolute scale for the T>0 transition was calculated. Since the QCD transition is a non- 
singular cross-over there is no unique T^. This well-known phenomenon was illustrated on the 
water- vapor phase diagram (the broadness of a transition was illustrated by the example of butter). 
Different observables led to different numerical values in the continuum and thermodynamic 
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limit also in QCD. Three observables were used to determine the corresponding transition temper- 
atures. The peak of the renormalized chiral susceptibility predicted rc=151(3)(3) MeV, whereas 
Tc based on the strange quark number susceptibility resulted in 24(4) MeV larger value. Another 
quantity, which is related to the deconfining phase transition in the large quark mass limit is the 
Polyakov loop. Its behavior predicted a 25(4) MeV larger transition temperature, than that of the 
chiral susceptibility. Another consequence of the cross-over are the non-vanishing widths of the 
peaks even in the thermodynamic limit, which were also determined. For the chiral susceptibil- 
ity, strange quark number susceptibility and Polyakov-loop we obtained widths of 28(5)(1) MeV, 
42(4)(1) MeV and 38(5)(1) MeV, respectively. 

These features, numbers and functions are attempted to be the full result for the T^Q transition, 
though other lattice fermion formulations — e.g. Wilson fermions (for ongoing projects see e.g. 
[ plj ^ ) or chiral fermions (for an early dynamical overlap test see [23], for the domain wall 
approach a recent presentation can be found in Ref . [^4|]) — are needed to cross-check the findings 
with staggered fermions. 

Results for the equation of state on lattices with A/f =4 and 6 were presented. Clearly, one should 
carry out the calculations on lattices with smaller lattice spacings and approach the continuum limit. 

Results at non- vanishing chemical potentials are even further from the continuum limit. At one 
single lattice spacings we determined the critical endpoint and saw indications for a third phase 
on the r-/i plane at large baryonic densities. As we emphasized in the abstract of our critical 
endpoint paper [17| „the continuum extrapolation is still missing". Using the present methods and 
technologies one can determine the continuum curvature for the T-pL phase diagram at /x=0 in 
the next few years. Unfortunately, using the present techniques the CPU capacity will be most 
probably not enough to carry out a controlled continuum extrapolation for the critical endpoint and 
for a conjectured colour superconducting phase in the next few years. 
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